Metagenome diversity illuminates the origins of pathogen effectors.

Recent metagenome-assembled genome (MAG) analyses have profoundly impacted Rickettsiology systematics. The discovery of basal lineages (novel families Mitibacteraceae and Athabascaceae) with predicted extracellular lifestyles exposed an evolutionary timepoint for the transition to host dependency, which seemingly occurred independent of mitochondrial evolution. Notably, these basal rickettsiae carry the Rickettsiales vir homolog (rvh) type IV secretion system and purportedly use rvh to kill congener microbes rather than parasitize host cells as described for later-evolving rickettsial pathogens. MAG analysis also substantially increased diversity for the genus Rickettsia and delineated a sister lineage (the novel genus Tisiphia) that stands to inform on the emergence of human pathogens from protist and invertebrate endosymbionts. Herein, we probed Rickettsiales MAG and genomic diversity for the distribution of Rickettsia rvh effectors to ascertain their origins. A sparse distribution of most Rickettsia rvh effectors outside of Rickettsiaceae lineages illuminates unique rvh evolution from basal extracellular species and other rickettsial families. Remarkably, nearly every effector was found in multiple divergent forms with variable architectures, indicating profound roles for gene duplication and recombination in shaping effector repertoires in Rickettsia pathogens. Lateral gene transfer plays a prominent role in shaping the rvh effector landscape, as evinced by the discovery of many effectors on plasmids and conjugative transposons, as well as pervasive effector gene exchange between Rickettsia and Legionella species. Our study exemplifies how MAGs can yield insight into pathogen effector origins, particularly how effector architectures might become tailored to the discrete host cell functions of different eukaryotic hosts.IMPORTANCEWhile rickettsioses are deadly vector-borne human diseases, factors distinguishing Rickettsia pathogens from the innumerable bevy of environmental rickettsial endosymbionts remain lacking. Recent metagenome-assembled genome (MAG) studies revealed evolutionary timepoints for rickettsial transitions to host dependency. The rvh type IV secretion system was likely repurposed from congener killing in basal extracellular species to parasitizing host cells in later-evolving pathogens. Our analysis of MAG diversity for over two dozen rvh effectors unearthed their presence in some non-pathogens. However, most effectors were found in multiple divergent forms with variable architectures, indicating gene duplication and recombination-fashioned effector repertoires of Rickettsia pathogens. Lateral gene transfer substantially shaped pathogen effector arsenals, evinced by the discovery of effectors on plasmids and conjugative transposons, as well as pervasive effector gene exchanges between Rickettsia and Legionella species. Our study exemplifies how MAGs yield insight into pathogen effector origins and evolutionary processes tailoring effectors to eukaryotic host cell biology.

Structural determination of Rickettsia lipid A without chemical extraction confirms shorter acyl chains in later-evolving spotted fever group pathogens.

Rickettsiae are Gram-negative obligate intracellular parasites of numerous eukaryotes. Human pathogens of the transitional group (TRG), typhus group (TG), and spotted fever group (SFG) rickettsiae infect blood-feeding arthropods, have dissimilar clinical manifestations, and possess unique genomic and morphological attributes. Lacking glycolysis, rickettsiae pilfer numerous metabolites from the host cytosol to synthesize peptidoglycan and lipopolysaccharide (LPS). For LPS, O-antigen immunogenicity varies between SFG and TG pathogens; however, lipid A proinflammatory potential is unknown. We previously demonstrated that Rickettsia akari (TRG), Rickettsia typhi (TG), and Rickettsia montanensis (SFG) produce lipid A with long 2' secondary acyl chains (C16 or C18) compared to short 2' secondary acyl chains (C12) in Rickettsia rickettsii (SFG) lipid A. To further probe this structural heterogeneity and estimate a time point when shorter 2' secondary acyl chains originated, we generated lipid A structures for two additional SFG rickettsiae (Rickettsia rhipicephali and Rickettsia parkeri) utilizing fast lipid analysis technique adopted for use with tandem mass spectrometry (FLATn). FLATn allowed analysis of lipid A structure directly from host cell-purified bacteria, providing a substantial improvement over lipid A chemical extraction. FLATn-derived structures indicate SFG rickettsiae diverging after R. rhipicephali evolved shorter 2' secondary acyl chains. While 2' secondary acyl chain lengths do not distinguish Rickettsia pathogens from non-pathogens, in silico analyses of Rickettsia LpxL late acyltransferases revealed discrete active sites and hydrocarbon rulers for long versus short 2' secondary acyl chain addition. Our collective data warrant determining Rickettsia lipid A inflammatory potential and how structural heterogeneity impacts lipid A-host receptor interactions.IMPORTANCEDeforestation, urbanization, and homelessness lead to spikes in Rickettsioses. Vector-borne human pathogens of transitional group (TRG), typhus group (TG), and spotted fever group (SFG) rickettsiae differ by clinical manifestations, immunopathology, genome composition, and morphology. We previously showed that lipid A (or endotoxin), the membrane anchor of Gram-negative bacterial lipopolysaccharide (LPS), structurally differs in Rickettsia rickettsii (later-evolving SFG) relative to Rickettsia montanensis (basal SFG), Rickettsia typhi (TG), and Rickettsia akari (TRG). As lipid A structure influences recognition potential in vertebrate LPS sensors, further assessment of Rickettsia lipid A structural heterogeneity is needed. Here, we sidestepped the difficulty of ex vivo lipid A chemical extraction by utilizing fast lipid analysis technique adopted for use with tandem mass spectrometry, a new procedure for generating lipid A structures directly from host cell-purified bacteria. These data confirm that later-evolving SFG pathogens synthesize structurally distinct lipid A. Our findings impact interpreting immune responses to different Rickettsia pathogens and utilizing lipid A adjuvant or anti-inflammatory properties in vaccinology.

Divergent Pseudomonas aeruginosa LpxO enzymes perform site-specific lipid A 2-hydroxylation.

Pseudomonas aeruginosa can survive in a myriad of environments, partially due to modifications of its lipid A, the membrane anchor of lipopolysaccharide. We previously demonstrated that divergent late acyltransferase paralogs, HtrB1 and HtrB2, add acyloxyacyl laurate to lipid A 2- and 2'-acyl chains, respectively. The genome of P. aeruginosa also has genes which encode two dioxygenase enzymes, LpxO1 and LpxO2, that individually hydroxylate a specific secondary laurate. LpxO1 acts on the 2'-acyloxyacyl laurate (added by HtrB2), whereas LpxO2 acts on the 2-acyloxyacyl laurate (added by HtrB1) in a site-specific manner. Furthermore, while both enzyme pairs are evolutionarily linked, phylogenomic analysis suggests the LpxO1/HtrB2 enzyme pair as being of ancestral origin, present throughout the Pseudomonas lineage, whereas the LpxO2/HtrB1 enzyme pair likely arose via horizontal gene transfer and has been retained in P. aeruginosa over time. Using a murine pulmonary infection model, we showed that both LpxO1 and LpxO2 enzymes are functional in vivo, as direct analysis of in vivo lipid A structure from bronchoalveolar lavage fluid revealed 2-hydroxylated lipid A. Gene expression analysis reveals increased lpxO2 but unchanged lpxO1 expression in vivo, suggesting differential regulation of these enzymes during infection. We also demonstrate that loss-of-function mutations arise in lpxO1 and lpxO2 during chronic lung infection in people with cystic fibrosis (CF), indicating a potential role for pathogenesis and airway adaptation. Collectively, our study characterizes lipid A 2-hydroxylation during P. aeruginosa airway infection that is regulated by two distinct lipid A dioxygenase enzymes.IMPORTANCEPseudomonas aeruginosa is an opportunistic pathogen that causes severe infection in hospitalized and chronically ill individuals. During infection, P. aeruginosa undergoes adaptive changes to evade host defenses and therapeutic interventions, increasing mortality and morbidity. Lipid A structural alteration is one such change that P. aeruginosa isolates undergo during chronic lung infection in CF. Investigating genetic drivers of this lipid A structural variation is crucial in understanding P. aeruginosa adaptation during infection. Here, we describe two lipid A dioxygenases with acyl-chain site specificity, each with different evolutionary origins. Further, we show that loss of function in these enzymes occurs in CF clinical isolates, suggesting a potential pathoadaptive phenotype. Studying these bacterial adaptations provides insight into selection pressures of the CF airway on P. aeruginosa phenotypes that persist during chronic infection. Understanding these adaptive changes may ultimately provide clinicians better control over bacterial populations during chronic infection.

Orientia tsutsugamushi: comprehensive analysis of the mobilome of a highly fragmented and repetitive genome reveals the capacity for ongoing lateral gene transfer in an obligate intracellular bacterium.

IMPORTANCE: Obligate intracellular bacteria, or those only capable of growth inside other living cells, have limited opportunities for horizontal gene transfer with other microbes due to their isolated replicative niche. The human pathogen Ot, an obligate intracellular bacterium causing scrub typhus, encodes an unusually high copy number of a ~40 gene mobile genetic element that typically facilitates genetic transfer across microbes. This proliferated element is heavily degraded in Ot and previously assumed to be inactive. Here, we conducted a detailed analysis of this element in eight Ot strains and discovered two strains with at least one intact copy. This implies that the element is still capable of moving across Ot populations and suggests that the genome of this bacterium may be even more dynamic than previously appreciated. Our work raises questions about intracellular microbial evolution and sounds an alarm for gene-based efforts focused on diagnosing and combatting scrub typhus.

Structural determination of Rickettsia lipid A without chemical extraction confirms shorter acyl chains in later-evolving Spotted Fever Group pathogens.

Rickettsiae are Gram-negative obligate intracellular parasites of numerous eukaryotes. Human pathogens of the Transitional Group (TRG), Typhus Group (TG), and Spotted Fever Group (SFG) rickettsiae infect blood-feeding arthropods, have dissimilar clinical manifestations, and possess unique genomic and morphological attributes. Lacking glycolysis, rickettsiae pilfer numerous metabolites from host cytosol to synthesize peptidoglycan and lipopolysaccharide (LPS). For LPS, O-antigen immunogenicity varies between SFG and TG pathogens; however, lipid A proinflammatory potential is unknown. We previously demonstrated that R. akari (TRG), R. typhi (TG), and R. montanensis (SFG) produce lipid A with long 2' secondary acyl chains (C16 or C18) compared to short 2' secondary acyl chains (C12) in R. rickettsii (SFG) lipid A. To further probe this structural heterogeneity and estimate a time point when shorter 2' secondary acyl chains originated, we generated lipid A structures for two additional SFG rickettsiae ( R. rhipicephali and R. parkeri ) utilizing Fast Lipid Analysis Technique adopted for use with tandem mass spectrometry (FLAT n ). FLAT n allowed analysis of lipid A structure directly from host cell-purified bacteria, providing substantial improvement over lipid A chemical extraction. FLAT n -derived structures indicate SFG rickettsiae diverging after R. rhipicephali evolved shorter 2' secondary acyl chains. Bioinformatics analysis of Rickettsia LpxL late acyltransferases revealed discrete active sites and hydrocarbon rulers for long versus short 2' secondary acyl chain addition. While the significance of different lipid A structures for diverse Rickettsia pathogens is unknown, our success using FLAT n will facilitate determining how structural heterogeneity impacts interactions with host lipid A receptors and overall inflammatory potential. IMPORTANCE: Deforestation, urbanization, and homelessness lead to spikes in Rickettsioses. Vector-borne human pathogens of Transitional Group (TRG), Typhus Group (TG), and Spotted Fever Group (SFG) rickettsiae differ by clinical manifestations, immunopathology, genome composition, and morphology. We previously showed that lipid A (or endotoxin), the membrane anchor of Gram-negative bacterial lipopolysaccharide (LPS), structurally differs in R. rickettsii (later-evolving SFG) relative to R. montanensis (basal SFG), R. typhi (TG), and R. akari (TRG). As lipid A structure influences recognition potential in vertebrate LPS sensors, further assessment of Rickettsia lipid A structural heterogeneity is needed. Here, we sidestepped the difficulty of ex vivo lipid A chemical extraction by utilizing FLAT n , a new procedure for generating lipid A structures directly from host cell-purified bacteria. These data confirm later-evolving SFG pathogens synthesize structurally distinct lipid A. Our findings impact interpreting immune responses to different Rickettsia pathogens and utilizing lipid A adjuvant or anti-inflammatory properties in vaccinology.

Orientia tsutsugamushi: analysis of the mobilome of a highly fragmented and repetitive genome reveals ongoing lateral gene transfer in an obligate intracellular bacterium.

The rickettsial human pathogen Orientia tsutsugamushi (Ot) is an obligate intracellular Gram-negative bacterium with one of the most highly fragmented and repetitive genomes of any organism. Around 50% of its ~2.3 Mb genome is comprised of repetitive DNA that is derived from the highly proliferated Rickettsiales amplified genetic element (RAGE). RAGE is an integrative and conjugative element (ICE) that is present in a single Ot genome in up to 92 copies, most of which are partially or heavily degraded. In this report, we analysed RAGEs in eight fully sequenced Ot genomes and manually curated and reannotated all RAGE-associated genes, including those encoding DNA mobilisation proteins, P-type (vir) and F-type (tra) type IV secretion system (T4SS) components, Ankyrin repeat- and tetratricopeptide repeat-containing effectors, and other piggybacking cargo. Originally, the heavily degraded Ot RAGEs led to speculation that they are remnants of historical ICEs that are no longer active. Our analysis, however, identified two Ot genomes harbouring one or more intact RAGEs with complete F-T4SS genes essential for mediating ICE DNA transfer. As similar ICEs have been identified in unrelated rickettsial species, we assert that RAGEs play an ongoing role in lateral gene transfer within the Rickettsiales. Remarkably, we also identified in several Ot genomes remnants of prophages with no similarity to other rickettsial prophages. Together these findings indicate that, despite their obligate intracellular lifestyle and host range restricted to mites, rodents and humans, Ot genomes are highly dynamic and shaped through ongoing invasions by mobile genetic elements and viruses.

Orientia and Rickettsia: different flowers from the same garden.

Recent discoveries of basal extracellular Rickettsiales have illuminated divergent evolutionary paths to host dependency in later-evolving lineages. Family Rickettsiaceae, primarily comprised of numerous protist- and invertebrate-associated species, also includes human pathogens from two genera, Orientia and Rickettsia. Once considered sister taxa, these bacteria form distinct lineages with newly appreciated lifestyles and morphological traits. Contrasting other rickettsial human pathogens in Family Anaplasmataceae, Orientia and Rickettsia species do not reside in host-derived vacuoles and lack glycolytic potential. With only a few described mechanisms, strategies for commandeering host glycolysis to support cytosolic growth remain to be discovered. While regulatory systems for this unique mode of intracellular parasitism are unclear, conjugative transposons unique to Orientia and Rickettsia species provide insights that are critical for determining how these obligate intracellular pathogens overtake eukaryotic cytosol.

Metagenome diversity illuminates origins of pathogen effectors.

Recent metagenome assembled genome (MAG) analyses have profoundly impacted Rickettsiology systematics. Discovery of basal lineages (Mitibacteraceae and Athabascaceae) with predicted extracellular lifestyles reveals an evolutionary timepoint for the transition to host dependency, which occurred independent of mitochondrial evolution. Notably, these basal rickettsiae carry the Rickettsiales vir homolog (rvh) type IV secretion system (T4SS) and purportedly use rvh to kill congener microbes rather than parasitize host cells as described for derived rickettsial pathogens. MAG analysis also substantially increased diversity for genus Rickettsia and delineated a basal lineage (Tisiphia) that stands to inform on the rise of human pathogens from protist and invertebrate endosymbionts. Herein, we probed Rickettsiales MAG and genomic diversity for the distribution of Rickettsia rvh effectors to ascertain their origins. A sparse distribution of most Rickettsia rvh effectors outside of Rickettsiaceae lineages indicates unique rvh evolution from basal extracellular species and other rickettsial families. Remarkably, nearly every effector was found in multiple divergent forms with variable architectures, illuminating profound roles for gene duplication and recombination in shaping effector repertoires in Rickettsia pathogens. Lateral gene transfer plays a prominent role shaping the rvh effector landscape, as evinced by the discover of many effectors on plasmids and conjugative transposons, as well as pervasive effector gene exchange between Rickettsia and Legionella species. Our study exemplifies how MAGs can provide incredible insight on the origins of pathogen effectors and how their architectural modifications become tailored to eukaryotic host cell biology.

Cryptic Genes for Interbacterial Antagonism Distinguish Rickettsia Species Infecting Blacklegged Ticks From Other Rickettsia Pathogens.

Background: The genus Rickettsia (Alphaproteobacteria: Rickettsiales) encompasses numerous obligate intracellular species with predominantly ciliate and arthropod hosts. Notable species are pathogens transmitted to mammals by blood-feeding arthropods. Mammalian pathogenicity evolved from basal, non-pathogenic host-associations; however, some non-pathogens are closely related to pathogens. One such species, Rickettsia buchneri, is prevalent in the blacklegged tick, Ixodes scapularis. While I. scapularis transmits several pathogens to humans, it does not transmit Rickettsia pathogens. We hypothesize that R. buchneri established a mutualism with I. scapularis, blocking tick superinfection with Rickettsia pathogens. Methods: To improve estimates for assessing R. buchneri infection frequency in blacklegged tick populations, we used comparative genomics to identify an R. buchneri gene (REIS_1424) not present in other Rickettsia species present throughout the I. scapularis geographic range. Bioinformatic and phylogenomics approaches were employed to propose a function for the hypothetical protein (263 aa) encoded by REIS_1424. Results: REIS_1424 has few analogs in other Rickettsiales genomes and greatest similarity to non-Proteobacteria proteins. This cohort of proteins varies greatly in size and domain composition, possessing characteristics of Recombination hotspot (Rhs) and contact dependent growth inhibition (CDI) toxins, with similarity limited to proximal C-termini (~145 aa). This domain was named CDI-like/Rhs-like C-terminal toxin (CRCT). As such proteins are often found as toxin-antidote (TA) modules, we interrogated REIS_1423 (151 aa) as a putative antidote. Indeed, REIS_1423 is similar to proteins encoded upstream of CRCT domain-containing proteins. Accordingly, we named these proteins CDI-like/Rhs-like C-terminal toxin antidotes (CRCA). R. buchneri expressed both REIS_1423 and REIS_1424 in tick cell culture, and PCR assays showed specificity for R. buchneri over other rickettsiae and utility for positive detection in three tick populations. Finally, phylogenomics analyses uncovered divergent CRCT/CRCA modules in varying states of conservation; however, only R. buchneri and related Tamurae/Ixodes Group rickettsiae carry complete TA modules. Conclusion: We hypothesize that Rickettsia CRCT/CRCA modules circulate in the Rickettsia mobile gene pool, arming rickettsiae for battle over arthropod colonization. While its functional significance remains to be tested, R. buchneri CRCT/CRCA serves as a marker to positively identify infection and begin deciphering the role this endosymbiont plays in the biology of the blacklegged tick.

Genome sequencing of Pseudomonas aeruginosa strain M2 illuminates traits of an opportunistic pathogen of burn wounds.

Pseudomonas aeruginosa is a Gram-negative nosocomial pathogen and one of the most prevalent organisms isolated from burn wounds worldwide. Pseudomonas aeruginosa strain M2 (O5 serotype, type B flagella) is utilized for examining the murine model associated with burns. Pseudomonas aeruginosa M2 is similar in lethality to common laboratory P. aeruginosa strains when infecting CD-1 mice. Conversely, we recently showed that, relative to these strains, P. aeruginosa M2-infected mice are more susceptible to sepsis and demonstrate a 6-log reduction in LD50 from subcutaneous infection at the infection site directly after 10% total body surface area burn. To better understand this striking phenotypic difference from other P. aeruginosa strains employed in burn models, we sequenced the P. aeruginosa M2 genome. A total of 4,136,641 read pairs were obtained, providing an average genome coverage of 97.5X; subsequent assembly yielded a draft genome with 187 contigs comprising 6,360,304 bp with a G + C content of 66.45%. Genome-based phylogeny estimation of 92 P. aeruginosa strains placed P. aeruginosa M2 with P. aeruginosa-12-4-4(59), a nonairway clinical strain isolated from the blood culture of a burn patient. Phylogenomic analyses identified genes shared between P. aeruginosa M2 and P. aeruginosa 14, another strain exhibiting increased lethality in thermal tissues, as well as P. aeruginosa M2 unique genes with diverse functions like degradation of toxic aromatic compounds, iron scavenging, swarming motility and biofilm formation, defense against invasive DNA, and host assault. Predicted lateral gene transfers illuminate proteins heretofore uncharacterized for roles in P. aeruginosa biology. Our work yields a rich resource for assessing P. aeruginosa genes required for increased lethality in burn tissue seroma.

Discovery of a Diverse Set of Bacteria That Build Their Cell Walls without the Canonical Peptidoglycan Polymerase aPBP.

Peptidoglycan (PG) is a highly cross-linked peptide-glycan mesh that confers structural rigidity and shape to most bacterial cells. Polymerization of new PG is usually achieved by the concerted activity of two membrane-bound machineries, class-A penicillin binding proteins (aPBPs) and class-B penicillin binding proteins (bPBPs) in complex with shape, elongation, division, and sporulation (SEDS) proteins. Here, we have identified four phylogenetically distinct groups of bacteria that lack any identifiable aPBPs. We performed experiments on a panel of species within one of these groups, the Rickettsiales, and found that bacteria lacking aPBPs build a PG-like cell wall with minimal abundance and rigidity relative to cell walls of aPBP-containing bacteria. This reduced cell wall may have evolved to minimize the activation of host responses to pathogens and endosymbionts while retaining the minimal PG-biosynthesis machinery required for cell elongation and division. We term these "peptidoglycan-intermediate" bacteria, a cohort of host-associated species that includes some human pathogens. IMPORTANCE Peptidoglycan (PG) is a large, cross-linked polymer that forms the cell wall of most bacterial species and confers shape, rigidity, and protection from osmotic shock. It is also a potent stimulator of the immune response in animals. PG is normally polymerized by two groups of enzymes, aPBPs and bPBPs working together with shape, elongation, division, and sporulation (SEDS) proteins. We have identified a diverse set of host-associated bacteria that have selectively lost aPBP genes while retaining bPBP/SEDS and show that some of these build a minimal PG-like structure. It is expected that these minimal cell walls built in the absence of aPBPs improve the evolutionary fitness of host-associated bacteria, potentially through evasion of PG-recognition by the host immune system.

Lipid A Structural Divergence in Rickettsia Pathogens.

Species of Rickettsia (Alphaproteobacteria: Rickettsiales) are obligate intracellular parasites of a wide range of eukaryotes, with recognized arthropod-borne human pathogens belonging to the transitional group (TRG), typhus group (TG), and spotted fever group (SFG) rickettsiae. Growing in the host cytosol, rickettsiae pilfer numerous metabolites to make a typical Gram-negative bacterial cell envelope. The O-antigen of rickettsial lipopolysaccharide (LPS) is immunogenic and has been shown to tether the S-layer to the rickettsial surface; however, little is known about the structure and immunogenicity of the Rickettsia lipid A moiety. The structure of lipid A, the membrane anchor of LPS, affects the ability of this molecule to interact with components of the host innate immune system, specifically the MD-2/TLR4 receptor complex. To dissect the host responses that can occur during Rickettsia in vitro and in vivo infection, structural analysis of Rickettsia lipid A is needed. Lipid A was extracted from four Rickettsia species and structurally analyzed. R. akari (TRG), R. typhi (TG), and R. montanensis (SFG) produced a similar structure, whereas R. rickettsii (SFG) altered the length of a secondary acyl group. While all structures have longer acyl chains than known highly inflammatory hexa-acylated lipid A structures, the R. rickettsii modification should differentially alter interactions with the hydrophobic internal pocket in MD2. The significance of these characteristics toward inflammatory potential as well as membrane dynamics between arthropod and vertebrate cellular environments warrants further investigation. Our work adds lipid A to the secretome and O-antigen as variable factors possibly correlating with phenotypically diverse rickettsioses.IMPORTANCE Spikes in rickettsioses occur as deforestation, urbanization, and homelessness increase human exposure to blood-feeding arthropods. Still, effective Rickettsia vaccines remain elusive. Recent studies have determined that Rickettsia lipopolysaccharide anchors the protective S-layer to the bacterial surface and elicits bactericidal antibodies. Furthermore, growing immunological evidence suggests vertebrate sensors (MD-2/TLR4 and noncanonical inflammasome) typically triggered by the lipid A portion of lipopolysaccharide are activated during Rickettsia infection. However, the immunopotency of Rickettsia lipid A is unknown due to poor appreciation for its structure. We determined lipid A structures for four distinct rickettsiae, revealing longer acyl chains relative to highly inflammatory bacterial lipid A. Surprisingly, lipid A of the Rocky Mountain spotted fever agent deviates in structure from other rickettsiae. Thus, lipid A divergence may contribute to variable disease phenotypes, sounding an alarm for determining its immunopotency and possible utility (i.e., as an adjuvant or anti-inflammatory) for development of more prudent rickettsiacidal therapies.

Evolution of Wolbachia mutualism and reproductive parasitism: insight from two novel strains that co-infect cat fleas.

Wolbachiae are obligate intracellular bacteria that infect arthropods and certain nematodes. Usually maternally inherited, they may provision nutrients to (mutualism) or alter sexual biology of (reproductive parasitism) their invertebrate hosts. We report the assembly of closed genomes for two novel wolbachiae, wCfeT and wCfeJ, found co-infecting cat fleas (Ctenocephalides felis) of the Elward Laboratory colony (Soquel, CA, USA). wCfeT is basal to nearly all described Wolbachia supergroups, while wCfeJ is related to supergroups C, D and F. Both genomes contain laterally transferred genes that inform on the evolution of Wolbachia host associations. wCfeT carries the Biotin synthesis Operon of Obligate intracellular Microbes (BOOM); our analyses reveal five independent acquisitions of BOOM across the Wolbachia tree, indicating parallel evolution towards mutualism. Alternately, wCfeJ harbors a toxin-antidote operon analogous to the wPip cinAB operon recently characterized as an inducer of cytoplasmic incompatibility (CI) in flies. wCfeJ cinB and three adjacent genes are collectively similar to large modular toxins encoded in CI-like operons of certain Wolbachia strains and Rickettsia species, signifying that CI toxins streamline by fission of large modular toxins. Remarkably, the C. felis genome itself contains two CI-like antidote genes, divergent from wCfeJ cinA, revealing episodic reproductive parasitism in cat fleas and evidencing mobility of CI loci independent of WO-phage. Additional screening revealed predominant co-infection (wCfeT/wCfeJ) amongst C. felis colonies, though fleas in wild populations mostly harbor wCfeT alone. Collectively, genomes of wCfeT, wCfeJ, and their cat flea host supply instances of lateral gene transfers that could drive transitions between parasitism and mutualism.

Risk1, a Phosphatidylinositol 3-Kinase Effector, Promotes Rickettsia typhi Intracellular Survival.

To establish a habitable intracellular niche, various pathogenic bacteria secrete effectors that target intracellular trafficking and modulate phosphoinositide (PI) metabolism. Murine typhus, caused by the obligate intracellular bacterium Rickettsia typhi, remains a severe disease in humans. However, the mechanisms by which R. typhi effector molecules contribute to internalization by induced phagocytosis and subsequent phagosomal escape into the cytosol to facilitate the intracellular growth of the bacteria remain ill-defined. Here, we characterize a new molecule, Risk1, as a phosphatidylinositol 3-kinase (PI3K) secreted effector and the first bacterial secretory kinase with both class I and III PI3K activities. Inactivation of Risk1 PI3K activities reduced the phosphorylation of phosphatidylinositol 4,5-bisphosphate to phosphatidylinositol 3,4,5-trisphosphate within the host, which consequently diminished host colonization by R. typhi During infection, Risk1 targets the Rab5-EEA1-phosphatidylinositol 3-phosphate [PI(3)P] signaling axis to promote bacterial phagosomal escape. Subsequently, R. typhi undergoes ubiquitination and induces host autophagy; however, maturation to autolysosomes is subverted to support intracellular growth. Intriguingly, only enzymatically active Risk1 binds the Beclin-1 core complex and contributes to R. typhi-induced autophagosome formation. In sum, our data suggest that Risk1, with dual class I and class III PI3K activities, alters host PI metabolism and consequently subverts intracellular trafficking to facilitate intracellular growth of R. typhiIMPORTANCERickettsia species are Gram-negative obligate intracellular bacteria that infect a wide range of eukaryotes and vertebrates. In particular, human body louse-borne Rickettsia prowazekii and flea-borne Rickettsia typhi have historically plagued humankind and continue to reemerge globally. The unavailability of vaccines and limited effectiveness of antibiotics late in infection place lethality rates up to 30%, highlighting the need to elucidate the mechanisms of Rickettsia pathogenicity in greater detail. Here, we characterize a new effector, Risk1, as a secreted phosphatidylinositol 3-kinase (PI3K) with unique dual class I and class III activities. Risk1 is required for host colonization, and its vacuolar phosphatidylinositol 3-phosphate generation modulates endosomal trafficking to arrest autophagosomal maturation. Collectively, Risk1 facilitates R. typhi growth by altering phosphoinositide metabolism and subverting intracellular trafficking.

A chromosome-level assembly of the cat flea genome uncovers rampant gene duplication and genome size plasticity.

BACKGROUND: Fleas (Insecta: Siphonaptera) are small flightless parasites of birds and mammals; their blood-feeding can transmit many serious pathogens (i.e., the etiological agents of bubonic plague, endemic and murine typhus). The lack of flea genome assemblies has hindered research, especially comparisons to other disease vectors. Accordingly, we sequenced the genome of the cat flea, Ctenocephalides felis, an insect with substantial human health and veterinary importance across the globe. RESULTS: By combining Illumina and PacBio sequencing of DNA derived from multiple inbred female fleas with Hi-C scaffolding techniques, we generated a chromosome-level genome assembly for C. felis. Unexpectedly, our assembly revealed extensive gene duplication across the entire genome, exemplified by ~ 38% of protein-coding genes with two or more copies and over 4000 tRNA genes. A broad range of genome size determinations (433-551 Mb) for individual fleas sampled across different populations supports the widespread presence of fluctuating copy number variation (CNV) in C. felis. Similarly, broad genome sizes were also calculated for individuals of Xenopsylla cheopis (Oriental rat flea), indicating that this remarkable "genome-in-flux" phenomenon could be a siphonapteran-wide trait. Finally, from the C. felis sequence reads, we also generated closed genomes for two novel strains of Wolbachia, one parasitic and one symbiotic, found to co-infect individual fleas. CONCLUSION: Rampant CNV in C. felis has dire implications for gene-targeting pest control measures and stands to complicate standard normalization procedures utilized in comparative transcriptomics analysis. Coupled with co-infection by novel Wolbachia endosymbionts-potential tools for blocking pathogen transmission-these oddities highlight a unique and underappreciated disease vector.

The complete mitochondrial genome of the cat flea, Ctenocephalides felis.

The cat flea, Ctenocephalides fells, is widely recognized as a global veterinary pest and a vector of pathogenic bacteria. We recently reported on the C. felis nuclear genome, which is characterized by over 38% protein coding gene duplication, extensive tRNA gene family expansion, and remarkable gene copy number variation (CNV) between individual fleas. Herein, we describe the assembly of the C. felis mitochondrial genome, a novel resource for comparative genomics of fleas and other insects. The order and content of mitochondrial genes is highly consistent with four previously sequenced flea mitochondrial genomes, limiting CNV to siphonapteran nuclear genomes.

Conjugative Transposons and Their Cargo Genes Vary across Natural Populations of Rickettsia buchneri Infecting the Tick Ixodes scapularis.

Rickettsia buchneri (formerly Rickettsia endosymbiont of Ixodes scapularis, or REIS) is an obligate intracellular endoparasite of the black-legged tick, the primary vector of Lyme disease in North America. It is noteworthy among the rickettsiae for its relatively large genome (1.8 Mb) and extraordinary proliferation of mobile genetic elements (MGEs), which comprise nearly 35% of its genome. Previous analysis of the R. buchneri genome identified several integrative conjugative elements named Rickettsiales amplified genomic elements (RAGEs); the composition of these RAGEs suggests that continued genomic invasions by MGEs facilitated the proliferation of rickettsial genes related to an intracellular lifestyle. In this study, we compare the genomic diversity at RAGE loci among sequenced rickettsiae that infect three related Ixodes spp., including two strains of R. buchneri and Rickettsia endosymbiont of Ixodes pacificus strain Humboldt, as well as a closely related species R. tamurae infecting Amblyomma testudinarium ticks. We further develop a novel multiplex droplet digital PCR assay and use it to quantify copy number ratios of chromosomal R. buchneri RAGE-A and RAGE-B to the single-copy gene gltA within natural populations of I. scapularis. Our results reveal substantial diversity among R. buchneri at these loci, both within individual ticks as well as in the I. scapularis population at large, demonstrating that genomic rearrangement of MGEs is an active process in these intracellular bacteria.

Rickettsia Lipid A Biosynthesis Utilizes the Late Acyltransferase LpxJ for Secondary Fatty Acid Addition.

Members of the Rickettsia genus are obligate intracellular, Gram-negative coccobacilli that infect mammalian and arthropod hosts. Several rickettsial species are human pathogens and are transmitted by blood-feeding arthropods. In Gram-negative parasites, the outer membrane (OM) sits at the nexus of the host-pathogen interaction and is rich in lipopolysaccharide (LPS). The lipid A component of LPS anchors the molecule to the bacterial surface and is an endotoxic agonist of Toll-like receptor 4 (TLR4). Despite the apparent importance of lipid A in maintaining OM integrity, as well as its inflammatory potential during infection, this molecule is poorly characterized in Rickettsia pathogens. In this work, we have identified and characterized new members of the recently discovered LpxJ family of lipid A acyltransferases in both Rickettsia typhi and Rickettsia rickettsii, the etiological agents of murine typhus and Rocky Mountain spotted fever, respectively. Our results demonstrate that these enzymes catalyze the addition of a secondary acyl chain (C14/C16) to the 3'-linked primary acyl chain of the lipid A moiety in the final steps of the Raetz pathway of lipid A biosynthesis. Since lipid A architecture is fundamental to bacterial OM integrity, we believe that rickettsial LpxJ may be important in maintaining membrane dynamics to facilitate molecular interactions at the host-pathogen interface that are required for adhesion and invasion of mammalian cells. This work contributes to our understanding of rickettsial outer membrane physiology and sets a foundation for further exploration of the envelope and its role in pathogenesis.IMPORTANCE Lipopolysaccharide (LPS) triggers an inflammatory response through the TLR4-MD2 receptor complex and inflammatory caspases, a process mediated by the lipid A moiety of LPS. Species of Rickettsia directly engage both extracellular and intracellular immunosurveillance, yet little is known about rickettsial lipid A. Here, we demonstrate that the alternative lipid A acyltransferase, LpxJ, from Rickettsia typhi and R. rickettsii catalyzes the addition of C16 fatty acid chains into the lipid A 3'-linked primary acyl chain, accounting for major structural differences relative to the highly inflammatory lipid A of Escherichia coli.

A Tangled Web: Origins of Reproductive Parasitism.

While typically a flea parasite and opportunistic human pathogen, the presence of Rickettsia felis (strain LSU-Lb) in the non-blood-feeding, parthenogenetically reproducing booklouse, Liposcelis bostrychophila, provides a system to ascertain factors governing not only host transitions but also obligate reproductive parasitism (RP). Analysis of plasmid pLbAR, unique to R. felis str. LSU-Lb, revealed a toxin-antitoxin module with similar features to prophage-encoded toxin-antitoxin modules utilized by parasitic Wolbachia strains to induce another form of RP, cytoplasmic incompatibility, in their arthropod hosts. Curiously, multiple deubiquitinase and nuclease domains of the large (3,841 aa) pLbAR toxin, as well the entire antitoxin, facilitated the detection of an assortment of related proteins from diverse intracellular bacteria, including other reproductive parasites. Our description of these remarkable components of the intracellular mobilome, including their presence in certain arthropod genomes, lends insight on the evolution of RP, while invigorating research on parasite-mediated biocontrol of arthropod-borne viral and bacterial pathogens.